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Abstract. It has been remarked that a fair measure of the impact of Atle 
Selberg's work is the number of mathematical terms which bear his name. One 
of these is the Selberg integral, an n-dimensional generalization of the Euler 
beta integral. We trace its sudden rise to prominence, initiated by a question 
to Selberg from Enrico Bombieri, more than thirty years after publication. In 
quick succession the Selberg integral was used to prove an outstanding con- 
jecture in random matrix theory, and cases of the Macdonald conjectures. It 
further initiated the study of q-analogues, which in turn enriched the Mac- 
donald conjectures. We review these developments and proceed to exhibit 
the sustained prominence of the Selberg integral, evidenced by its central role 
in random matrix theory, Calogero-Sutherland quantum many body systems, 
Knizhnik-Zamolodchikov equations, and multivariable orthogonal polynomial 
theory. 



1. Discovery and reappearance 

1941 and 1944. With the passing of Atle Selberg on August 6th 2007 at age 90, 
it is timely to reflect on his mathematical legacy. Indeed a number of brief articles 
highlighting some of his most influential mathematical discoveries were written 
shortly after the news of his death, see e.g., [74]. It is our aim to add to these 
tributes by giving a more comprehensive account of the mathematics, both pure 
and applied, related to what now is referred to as the Selberg integral 



(1.1) 5„(a,/3,7):= 




^ yr T{a + n)TiP + J7)r(l + (j + 1)7) 

The evaluation of this integral is valid for complex parameters a, /3, 7 such that 

(1.2) Re(a) > 0, Re(/?) > 0, Re(7) > - min{l/n, Re(a)/(n - 1), Re(/3)/(n - 1)}, 

corresponding to the domain of convergence of the integral. 

The proof of ( |1.1[ ) is the subject of Selberg's 1944 paper "Bemerkninger om et 
multipelt integral" (Remarks on a multiple integral) [134] — the only one of Sel- 
berg's works written in Norwegian — published in Norsk Matematisk Tidsskrift. 
The latter has been compared [22] to the Scandinavian equivalent of the Mathe- 
matical Gazette, with contents ranging from short research papers on subjects of 
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general interest to discussions on teaching problems. Selberg himself remarks in his 
collected works [136] that 

This paper was published with some hesitation, and in Norwegian, 
since I was rather doubtful that the results were new. The journal 
is one which is read by mathematics-teachers in the gymnasium, 
and the proof was written out in some detail so it should be under- 
standable to someone who knew a little about analytic functions 
and analytic continuation. 

Selberg's proof of ( |1.1[ ) proceeds by supposing 7 is a positive integer, and ex- 
panding 

l<i<j<n 0<fci,...,fe„<2(n-l)7 

Substituting this expansion in the definition of Sn{a, allows the resulting in- 
tegrals to be evaluated by the Euler beta integral [44] 



(1.3) 



Jo 



r(a)r(/j) 



which itself is (1.1 1 with n = 1. The details, in English, of the proof from here 



on can be found in [54,107], for example. Perhaps the most significant feature is 



the final step. It requires analytically continuing off the integers. Thus with (1.1) 
established for 7 a positive integer, the remaining task is to establish its validity 
for all complex 7 such that both sides are well defined. 



For this purpose, after noting that both the left- and right-hand side of ( 1.1 1 are 
bounded analytic functions of 7 for Re(7) > 1 at least, Carlson's theorem can be 
used [26]. The latter applies to functions f{z) analytic for Ke{z) > satisfying the 
bound \f{z)\ — O(e^l^l), fi < tt. The theorem asserts that if furthermore f{z) = 
on the nonnegative integers then, identically, f{z) = 0. Note that the function 
f{z) — sinTTz shows that the bound /x < tt is optimal. Selberg did not make direct 
use of Carlson's theorem, but rather derived from first principles the same result in 
the case that f{z) is bounded in the right half-plane, which is all that is required 



to finalise the proof of ( 1 . 1 1 



Interestingly, although [134] contains the first proof of (1.1 1, it is not the first 
time it appeared in print. This occurred three years earlier (albeit with the change 
of variables ti = Si/{1 + Si) so that Si S [0,C!o)) in Selberg's 1941 paper "Uber 
einen Satz von A. Gelfond" (On a theorem of A. Gelfond) [133]. Like [134], this 
earlier paper appeared in a Norwegian journal, this time Archiv for Mathematik og 
Naturvidenskab, known for having Sophus Lie as one of its founders. In a footnote 
Selberg remarks 

Leider habe ich die Formcl (11) [The Selberg integral] nirgends in 
der Litteratur finden konnen, ein Beweis hier zu bringen schcint 
aber nicht angebracht, da die Arbeit sonst zu sehr anschwellen 
wiirde; soUte sich aber herausstellen, dass die Formcl ncu ware, 
beabsichtige ich spater ein Beweis zu veroffentlichen. 

[Unfortunately I have been unable to find formula (11) [The Selberg 
integral] in the literature. To present a proof here, however, seems 
inappropriate, as it would make this paper significantly longer. If 
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it turns out that the formula is new, I intend to pubUsh a proof at 
a later date.] 

Curiously, Selberg used his integral in [133] to prove a result of some similarity to 
Carlson's theorem. As already noted, the latter is itself an ingredient in Selberg's 



proof of (1.1 1. Selberg's result relates to entire functions f{z) of exponential type 
cr(/), defined by 

cr(/) :=limsup-log( max|/(z)|). 

r^oo r V \z\=r J 

A theorem of Hardy and Polya [21,122] states that if cr(/) < log 2 and / takes integer 
values at the nonnegative integers, then /(z) is polynomial. The transcendental 
function /(z) = V shows that this bound is optimal. A. Gelfond [64] generalized 
this by proving that if (t(/) < nlog(l + exp(l/n — 1)) and /, together with its first 
n — 1 derivatives, take integer values at the nonnegative integers, then /(z) is a 
polynomial. However, for n > 1 this bound is not optimal. By using his integral 
Selberg improved Gelfond's bound for n > 1 to (t(/) < logm„, where to„ is the 
minimum value of Il"^]^(l + yi) under the conditions yi > 0, yi ■ ■ ■ yn ~ e^~" and 
ni<i<j<rJyr^ ~ ^7^1 ~ ^' "^^^^ improves Gelfond's result since 

n 

n(i + y.)>(i + e^/"-^)". 

i=l 

The 1950s to the late 1970s — the Mehta integral. For over thirty years the 
Selberg integral went essentially unnoticed. It was used only once — in the special 
case a = /3=l, 7 = 2 — in a study by S. Karlin and L.S. Shapley relating to the 
volume of a certain moment space, published in 1953 [87]. 



Around ten years later there was again good reason to make use of ( 1.1 ). Build- 
ing upon the earlier work of E.P. Wigner in the 1950s, F.J. Dyson wrote a series of 
papers on the statistical theory of energy levels of complex systems. These papers 
ranged from the theory's foundations to its practical use in the analysis of experi- 
mental data. This last topic was addressed in part V of the series, in a work, written 
jointly with M.L. Mehta and published in 1963, which also summarizes both the 
status of the theory and open problems from that date. 

A basic point is that random Hermitian matrices are used to model the highly 
excited states of complex nuclei. These matrices are taken to have real, complex 
or real quaternion elements, and correspond to the quantum system having time 
reversal symmetry, no time reversal symmetry, or time reversal symmetry with an 
odd number of spin 1 /2 particles respectively. 

In the real case all independent elements are chosen from independent standard 
normals N[0, 1], and in the complex case the diagonal elements are chosen indepen- 
dently from N[0, 1] while the off-diagonal elements are chosen independently from 
N[0, 1/V2] + i N[0, l/\/2]. Real quater nion elements are themselves 2x2 blocks of 
the form 

z 

—w 

In general, the eigenvalues of matrices with real quaternion elements are doubly 
degenerate. On the diagonal, each independent entry is chosen from N[0,-\/2], 
while on the off-diagonal, each independent entry is chosen from N[0, 1] -I- i N[0, 1]. 
These ensembles of random matrices are referred to as the Gaussian orthogonal, 
unitary and symplectic ensembles (abbreviated GOE, CUE and GSE) respectively. 
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For each of the three Gaussian ensembles the joint probability density function 
(PDF) for the eigenvalues can be computed explicitly as [107] 

n 

Here /? = 1, 2, 4 for the GOE, GUE and GSE respectively, while F„ is the normal- 
ization 

(1-5) ^»(7):=7^/ •••/ n^"*-/' n It,-^,rdti...dt„, 



referred to as Mehta's integral. In [108] Mehta and Dyson evaluated F„(/3/2) for 
each of the the three special values of (3. Combining this with the evaluations for 
n = 2 and n — i for general (3 led them to conjecture that 



(1.6) F„(7) = n ^ 



-jl) 



r(i + 7) 



Assuming the validity of (1.6 1 for general 7, certain averages associated with 



( 1 .4 1 at the special random matrix couplings /3 = 1,2,4 are accessible. This becomes 
apparent by writing ( 1 .4 1 in the form 

The mean ([/) for a given /3 is now computed by taking the logarithmic derivative of 
the normalization _F'„(/3/2). A further differentiation with respect to (3 then yields 
the fluctuation (C/^) - (U)^. 

The form (1.7 1 highlights an analogy with the equilibrium statistical mechanics 
of a classical gas of n particles on the line, at inverse temperature f3, interacting 
via a repulsive logarithmic Coulomb potential and confined by a harmonic well. 
The quantity exp(— /?[/) is referred to as the Boltzmann factor. This interpretation 
plays a prominent role in Dyson's series of works. Indeed, the notation for the 
averages used above stems from the statistical physics literature (and corresponds 
to the mean energy and specific heat of the Coulomb gas) and may be substituted 
by the mean ^{U) and variance cr'^{U) respectively. 

It is not hard to see that the Selberg integral can be used to evaluate Mehta's 
integral thus proving the conjecture (1.6 1. By the change of variables 1— s- (1 — 
U/L)/2 in (|rT} 

(1.8) lim 2^'(2L)"+"("-i)'^5„(LV2,LV2,7) - F^{^). 

L — >oo 



Use of Stirling's formula to compute the same limit on the right-hand side of ( 1.1 ) 
then gives (1.6 1. However, in 1963 when Mehta and Dyson published their conjec- 
ture the Selberg integral was essentially unknown and so this method of proof was 
not available. 

The Mehta-Dyson conjecture received more prominence with its appearance in 
the first edition of Mehta's book Random Matrices and the Statistical Theory of 
Energy Levels, published in 1967 [105]. In 1974 Mehta submitted the conjecture 
to the problems section of SIAM Review [106], thus gaining exposure to an even 
wider mathematical audience. A proof, exactly the one mentioned in the previous 
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paragraphs, was finally uncovered in the late 1970s by Selberg's IAS colleague 
Enrico Bombieri. The remarkable story behind this proof is best told in Bombieri's 
own words [23]: 

Since 1976 I had been studying elementary methods in prime num- 
ber theory and in particular a several variable extension of Cheby- 
shcv's well-known method to obtain upper and lower bounds for 
the number of primes up to a given bound. In the course of my re- 
searches I came across the problem of the asymptotic computation 
of certain multiple integrals, the simplest being 



where p, q, r are large positive integers and the integral is extended 

to the product of the unit circles |z,;| = 1 or to [0, 1]". 

The integral is related to a partition function for the one-di- 
mensional Coulomb gas on the unit circle \z\ = 1 with a fixed 
point charge at z = 1, as it was explained to me by my friend and 
colleague Tom Spencer, so I went to Dyson and asked him whether 
physicists had encountered such things before; maybe he could save 
me some efforts. 

Dyson told me that for q = 1 and 2 an integral of this type, 
over the real line with a gaussian measure, had indeed been studied 
and he referred me to a book by Mehta. Then I went to see Atle 
to ask his opinion about what I was doing in order to study the 
distribution of primes and whether he felt it was of any interest 
and whether he had any opinion on it. 

He immediately recognized my integral as a complex version of 
the generalized beta integral he had studied before and he gave me 
an off-print of his paper. It was not difficult to follow his proof, 
given for an integral over [0,1]", and use a classical method to 
write a Beta integral as a complex integral to solve my problem of 
computing my integral exactly. The multiple integral over [0. 1]" 
is of course Selberg's integral, as in that case arithmetical appli- 
cations require r to be a large negative (not positive) integer. It 
was also quite easy to get a confliient form of the Selbcrg integral 
and compute exactly the Mehta integrals for a general value of the 
parameter and make physicists happy. 

Since this was of interest to Dyson, I went back to Dyson and 
told him that using the Selberg integral one could compute the 
integral of interest to physicists. 

More from the 1960s and 70s — constant term identities. The consider- 
ation of time reversal symmetry leading to three ensembles of Hermitian matrices 
applies equally well to unitary matrices [42]. A conventional time reversal sym- 
metry requires that U = U^, no time reversal symmetry imposes no constraint, 
whilst a time reversal symmetry for a system with an odd number of spin 1 /2 par- 
ticles requires U = (here D denotes the quaternion dual; see e.g., [54, Ch. 2]). 
Choosing such matrices with a uniform probability then gives what are referred 




'n 
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to as the circular orthogonal, unitary and symplectic ensembles (COE, CUE and 
CSE) respectively. Their joint eigenvalue PDFs are given explicitly by 



(1.9) 



(27r)"C„(/?/2) 



where C„ is the normalization 

1 



(1.10) 



Cnil) = 



(27r)^ 



n Ie' 

1<? <j<n 



n I' 



-e'^^f^d0i---d0„, 



and /3 = 1,2,4 for the COE, CUE and CSE respectively. 



As for (1.5), the random matrix calculations give (1.10 1 in terms of gamma 



functions for the three special values of (i. Furthermore, the case n = 2 for general 



P can be related to the Euler beta integral (1.3 1, whilst the case n = 3 gives a sum 



which is a special instancec of an identity of Dixon for a well-poised 3F2 series [6], 



(1.11) 3F2 



a, b, c 
a - 6, 1 4 



r(i 



t)r(i 
a)r(i- 



a - b)T{l + a - c)r(l + I - 6 - c) 
f - b)T{l + f - c)r(l + a - 5 - c) ■ 



Based on all of these results, Dyson, in part I of his series of papers, made the 
conjecture [42] 



(1.12) 



Cnh) 



r(l + nj) 

r"(i + 7)' 



In the same paper, Dyson observed that with 7 a nonnegative integer, say k, (1.10) 



can be rewritten as the constant term (CT) in a Laurent expansion. This allows 

(fcn)! 



(1.12) to be rewritten as 



(1.13) 



OT n (-1;) (-S) 

l<i<j<n ■> 



{kiy 



This constant term identity, and thus, by Carlson's theorem, the conjecture (1.12) 



was soon proved by J. Gunson and K. Wilson [159], and later in a very efficient 
analysis by I.J. Good [65]. Gunson's proof is mentioned in [42], but the work is 
unpublished; reference often given to [66] in this context actually refers to the proof 
of another conjecture of Dyson. Twenty years after his proof Wilson was to receive 
the Nobel Prize in physics for his work on the renormalisation group approach to 
the study of critical phenomena; [159] is his first publication. 

In their proof, Wilson and Good both took advantage of the extra degrees of 
freedom offered by Dyson's more general conjecture, also contained in [42], 



(1.14) 



CT 



n(> 



(ai H h a„)! 



ai! 



The formulation of this was in turn motivated by the extra degrees of freedom 
permitted by Dixon's identity (1.11 1, to which (1.141 reduces in the case n = 3. 



In fact, as observed by R. Askey [11], the Selberg integral can be used to prove 
Dyson's conjecture (1.12) directly without the need for (1.14). Askey's observation 
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is based on the easily established general identity 

(1.15) / ■■■ [ {h---tr,)^-^f{h,...,tr,)dti---dU, 

Jo Jo 



V2sin<y 



,iC(ei+---+e„) /(-e'^i, 



-e'^")d0i---d0„, 



valid for / a Laurent polynomial and Re(C) large enough so that the left-hand side 
exists. Applying (1.15) to the Selberg integral with /? a positive integer and 7 a 

M„(a,6,7), 



nonnegative integer shows that 

(1.16) Sn{a,(3,j)^{-ir+i"^hf^y 

Vsm TTO/ 

b — [n— 1)7, (3 :— a + b+l and 



where a : 
(1.17) Af„(a,&,7): 



X n 

l<i<_7<n 



From (1.16 1, the Selberg integral, the reflection formula and finally Carlson's the- 
orem, it follows that 



(1.18) 



r(l + a + 6 + j7)r(l + (j + l)7) 



(l + a + j7)r(l + 6 + j7)r(l + 7)' 



for a, 6, 7 e C such that 

Y{_e{a + b+\) > 0, Re(7) > - min{l/n, Re(a -I- 6 -f l)/(n - 1)}. 



For a = b = this is Dyson's conjecture (1.12). 

The change of variables e'*'' = (i — ti)/{i +ti) maps the unit circle onto the real 
line via a stereographic projection. Applying this to the integral (1.17) leads to 



(1.19) 



(2 



•••/ TTt^ ^ TT K»-i,f^dii---dt 



2-ri(Q+/3-l)+ri(n-l)7 JJ^ 
J=0 



r(« + /? - 1 - (n + j - i)7)r(i + (.7 + 1)7) 

r(a-j7)r(/3-j7)r(l + 7) 



When n—1 this is the Cauchy beta integral. 

In the letter to Dyson reprinted on the next page Selberg communicated the 



multiple Cauchy integral ( 1.19 1. Subsequently, in a letter to Askey dated 25 March 
1980 , he mentioned both (1.17) and ( 1.19 1, and pointed out their exact relationship. 



The first time ( |1.17| appeared in print w as in W.G. Morris' 1982 PhD thesis [113]. 
In his thesis Morris provided a proof of (1.17) along the lines of Selberg's proof of 
(1.1 1, and applied it to obtain constant term identities. For these reasons ( 1.17[ ) is 
now commonly referred to as the Morris integral. 
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A culmination — The Macdonald Conjectures. In 1982 I.G. Macdonald [100] 



published his now famous ex-conjectures, generalizing the Mehta integral (1.5 1 to 



all finite reflection or Coxeter groups G, and the Dyson constant term identity 



(1.131 to all finite root systems. 

Let G be a finite group of isometrics of M" , generated by refiections in hyper- 
planes. Normalise (up to sign) so that the equations for the hyperplanes take the 
form 

anxi + • • • + ainXn — with a^^ + • • • + a^„ = 2, 
where i labels the hyperplanes, and form the polynomial 

N 

P{x) = W{a,iiXi H h ainXn)- 

i=l 

Geometrically, 2-^/'^P{x) gives the product of the distances of the point x S M to 
the N hyperplanes. 

By its action on M" the group G acts on polynomials in x = {xi, . . . ,Xn)- The 
polynomials that are invariant under the action of G are referred to as G-invariant 
polynomials. They form an M-algebra ...,/„] generated by n algebraically 

independent polynomials /i, . . . , /„ of degrees di, . . . , d„. Unhke the set of /^'s, the 
set of di's is uniquely determined by the underlying reflection group. 

A final ingredient required in the Macdonald integral conjectures is the Gaussian 
measure (p on M" 

e-l^l'/2 
■= (^^^da;i---dx„, 

where |a;p :— J27=i ^i- 

With the above notations Macdonald's (ex)-conjecture [100, Conjecture 5.1] 
states that for each finite refiection group G 



(1.20) y'|p(^)|27d^(:,)^ jjixi- 



(1 + 7) ■ 



For the three infinite families of crystallographic refiection groups (or refiection 
groups of Weyl type) A„_i, B„ and D„ the Macdonald conjecture follows as a limit 
of the Selberg integral. For type A„_i this corresponds to Bombieri's proof of the 
Mehta integral mentioned earlier. For types B„ and D„ this is due to A. Regev, 
although the actual proof appeared for the first time in the paper by Macdonald, 
to whom Regev communicated his results. 

Around the same time as Macdonald formulated his conjectures Regev was study- 
ing the large n behaviour of sums of the form Sg{n) :— X]a('^->^)^ where the sum is 
over partitions A of at most £ parts, and d\ is the dimension of the irreducible (3„ 
character [A] . Combining the hook- length formula for d\ with a careful asymptotic 
analysis, Regev showed [127] (see also [29]) that the asymptotics of sums like Sf{n) 
leads exactly to Mehta's integral. Regev remarks [128] 

From reactions to preprints and talks on [127], first from Richard 
Stanley (who in 1978 attended my seminar talk at UCSD) then 
from Freeman Dyson, I learned about the Mehta and the Macdon- 
ald conjectures. In a letter, Dyson also mentioned that the Mehta 
conjecture had just been verified — by applying the Selberg inte- 
gral. William Beckner then showed me the details of how to deduce 
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the Mehta — and some other integrals — from the Selberg integral. 
I worked on the other classical cases of the Macdonald conjecture 
and managed to verify these shortly afterwards, in 1979. 
The Coxeter group A„_i is the symmetry group of the (n — l)-simplex. It is a 
group of order n! generated by the n(n — l)/2 hyperplanes 

Xi — Xj ^ Q for ^ < i < j < n, 



and is isomorphic to the symmetric group ©„. All of the ingredients in (1.20) 
can thus easily be determined explicitly. The polynomial P{x) is given by the 
Vandermonde product 

(1.21) P{x)^ W [x^-Xj)^: ^{x) 

l<i<j<n 

and the G-invariant polynomials are given by the symmetric polynomials in x. 
One of the classical bases for the algebra of symmetric functions is given by the 
elementary symmetric functions {ei, . . . , e„} with 

6j.(x) — ^ ^ Xi^Xi2 ' ' ' Xi^. 



Accordingly the set of degrees di is given by {l,2,...,n}, and (1.20) reduces to 
Mehta's integral (1.5). 

The Coxeter groups B„ and D„ are the symmetry groups of the n-cube and 
n-demicube, and for these groups ( |1.20 1 takes the form 



l<i<j <n 



r(i + 7) 



and 



ni 2 2|27 n / X r(l + n7) 'i-r^r(l + 2i7) 

respectively. Making the changes ti — xf/ (2L), a = c + 1/2 and /3 = L + 1 in the 
Selberg integral and letting L tend to infinity gives 

g„ 1=1 i<i<j<n J=0 ^ > J U \ > U 

where on the left use has been made of Legendre's duplication formula. The above 
integral, known as the BC„ Mehta integral, leads to the B„ and D„ integrals by 
setting c — J and c = respectively. 

In his original paper Macdonald established several other instances of his conjec- 
ture, not relying on the Selberg integral. For 7=1 Macdonald presented a uniform 
proof for all crystallographic reflection groups. Another case of (1.20) — one that 
may be verified by purely elementary means — is that of the dihedral group l2(m), 
the symmetry group of a regular m-gon. 

A uniform proof of Macdonald's conjecture for all crystallographic reflection 
groups — A„_i, B„, D„, Eg, E7, Eg, F4 and G2 — was found in 1989 by E. Opdam 
[119] using the Heckman-Opdam theory of hypergeometric shift operators [71,119]. 
Several years later, combined theoretical and computer efforts by Opdam [120] and 
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F. Garvan [59, 60] also dealt with the remaining non-crystallographic groups H3 
and H4, thereby settling Macdonald's conjecture in full. 

In his paper Macdonald put forward many further conjectures related to root 
systems. One of these [100, Conjecture 2.7] is the generalization of Dyson's constant 



term identity (1.131 to arbitrary (finite) root systems. Let $ be a root system (not 
necessarily reduced) with corresponding Weyl group W. For a e $ let exp(a) be 
a formal exponential. Denote the degrees of the fundamental invariants of W by 
di, . . . ,di. The di may, for example, be obtained from the simple formula 

1 - t^tia)+s(a) ^ ' 1 _ ^d. 

11 l_iht(a) LL^—j-^ 
a6$+ 1=1 

where is the set of positive roots of the root system, ht(Q;) is the height of the 
root a and s{a) = 1 if a/2 ^ <i>+ and s{a) = 2 if a/2 e <i>^. (The latter can only 
occur for nonreduced root systems.) Then Macdonald's constant term conjecture 
asserts that 

(1.22) CT n(i-e")' =n 



k 



For the root system A„_i, $ = {e^ — ej|l < i, j < n, i ^ j} with the ith standard 
unit vector in M". The degrees in this case are given by 2,3, ... ,n so that, after 
the identification exp(ei — €j) = Xi/xj, one recovers Dyson's conjecture. 



When fc = 1 equation (1.221 simply follows from the classical Weyl denominator 
formula. Macdonald also proved the k = 2 case using algebraic techniques. Once 
again the Selberg integral implies the conjecture for all infinite series! 0^^, D^ 
and BC„. Since the first three are all contained in the latter the most succinct 
derivation arises by slightly generalizing the problem — Macdonald does this for 
all root systems in [100, Conjecture 2.3] — and considering the constant term of 

CT Yi (l-e")''". 

CtG'I'BCn 

Here 

$BC„ = {±q|1 <i<n}U {±2e,|l < i < n} U {±e, ±ej|l < i < j < n} 

is the BC„ root system and k^ = fci if a is of type zte^, k^ — k^ if a is of type 
±2ei and ka ~ k2 otherwise. The root systems B„, C„ and D„ are then obtained 
by taking fcs = or fci = or fci = fcs = respectively. By the substitution 
exp(ei) 1-^ exp(2i^?i) it follows that 

9n{a+6)+n(n— l)c 

(1.23) CT n (l-e")'=" = 

X / ••• / nsin»(0,)cos''(0,) H Sin^(0,-0,)sin^(0, + 0j)d0i-..d0„, 

■^0 -^0 »=1 l<t<j<n 

with a — 2ki + 2k^, b ~ 2k^ and c = 2^2. Introducing new integration variables 
ti — sin^(6'i) for all 1 < i < n the integral on the right transforms into the Selberg 
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integral, so that by (1.1 1 and the Legendre duphcation formula 
CT Y[ (l-e")'=° 

^ri(/£i+2fc3)+n(n — l)fe2 

Sn{ki + fcs + i, fcs + i, fc2) 



n 



(fc2 + ifc2)!(2fci + 2/c3 + 2ik2y.{2k3 + 2ik2)\ 



fj-^ fc2!(fci + fca + ik2)\{k3 + i/c2)!(A:i + 2^3 + (n + z - l)fc2)! 
In a not dissimilar manner D. Zeilberger [163] showed that the n — 2> case of the 



Morris integral (1.17) leads to the Macdonald conjecture for the exceptional root 
system G2. This result later found application in a study linking random matrix 
theory to number theoretical L-functions [91] (see also the section on the value 
distribution of log (,{1/2 + it) below). 

A unified proof of \\.22 ) for all root systems, based on hypergeometric shift 



operators, is again due to Opdam [119]. Pages 18 20 below contain an outline of 
this proof for the root system A„_i. 

2. Underpinnings of the Selberg integral 



The Dixon— Anderson integral. The Euler beta integral (1.3 1 has for its inte- 
grand the product of power functions x"'~^y^^^^ with y = 1 — x. It is evaluated 
as a ratio of gamma functions, which in turn are integrals over the product of a 
power function and exponential function. In the theory of finite fields, the role of 
power and exponential functions are played by multiplicative and additive charac- 
ters. These can be used to define the finite field analogue of the gamma and beta 
integrals, known as the Gauss and Jacobi sums respectively. Moreover, these finite 
field quantities satisfy an analogue of the beta integral. From Selberg's commen- 
tary [136], we know that in the 1940s he investigated finite field analogues of ([TT}, 
and formulated a conjecture which he could prove only for n = 2. The existence of 
such finite field analogues was revealed by Selberg in the letter to Askey dated 25 
March 1980, referred to on page|7] Selberg also mentioned this in some colloquium 
lectures. A member of the audience on one of these occasions, Ron Evans, has 
provided us with the following recollection [51]: 

Somewhere around 1980, Selberg came to UCSD for a colloquium 
talk. Some department members at the UCSD talk were shocked by 
the subject matter. They were expecting to hear about his recent 
work, but instead his entire talk was on the Selberg integral. I was 
fascinated to learn of this integral, and ended up writing several 
papers on q-analogues and on finite field analogues. One of these 
(published in 1981) formulated n-dimensional finite field analogues, 
which I was able to prove for n = 2. Selberg had mentioned in his 
talk that he had finite field analogues for n = 2, so I was reluctant 
to write up my proof. However, some people who knew Selberg 
told me that he'd never ever get around to publishing his proof, 
so I took the bold step of asking permission to include my proof 
for n = 2 with my general conjecture (with due credit, of course). 
He generously wrote back that he didn't mind if I publish a proof 
of the "right" version of the theorem, but that he didn't want to 
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be credited with my version, which was too weak! So I proved 
the stronger theorem for n — 2 that he supphed in his letter, and 
that led to stronger conjectures for general n (ultimately proved by 
Anderson [2]). 

The finite field paper of Evans referred to above is [46] , and Selberg in his commen- 
tary [136] references this as being his state of knowledge from the 1940s. In fact, 
the Anderson paper left open some of the conjectures from [46] and Evans himself 
was able to apply Anderson's approach to provide the remaining proofs [47]. For a 
detailed account of the finite field Selberg integral, we refer to [6] . 

In 1991, motivated by the quest for a proof of the finite field conjecture, G.W. An- 
derson [3] published a proof of the Selberg integral based on another multiple inte- 
gral, namely 

n n+1 

(2.24) / II {U-t,)l[l[\U-a,\'^~Uh---dtn 

•'^ l<i<j<n i=l j = l 



nrjir(^.) 

r(E 



n+1 



n 

l<i<j<n+l 



where X is the domain of integration ai > ti > a2 > t2 > ■ ■ ■ > tn > a„+i, and 
Re(si) > 1 for 1 < i < n -f 1. Anderson's idea was to use (2.241 to compute in two 
different ways the integral 



if(a,/3,7) := 



X' 



n+1 

n 

i=l 



71 n+1 

4 = 1 i = l 



17-1 



yj\ n 1^^ 

l<!<j<n+l 



X n \y'' 

l<z<j'<n 

where X' denotes the domain of integration 

1 > xi > yi > X2 > 2/2 > • ' 
First integrating over the y-variables gives 

r"+i(7) 



Xj\dxi--- da;„+id?/i • • • d?/„ , 



X(a,/3,7) 



> yn> > 0. 



S'„+i(a,/3,7), 



(n+ l)!r((n-f 1)7) 
while first integrating over the a;-variables gives 

K{a, (3, 7) = r(a)r(/3)r"(7) ^ y 

n\L(a + p + nj) 

Equating the two forms reveals a first order recurrence for the Selberg integral in 
n. Together with the initial condition So{a,P,j) = 1 this reclaims (1.1 1. 



A large portion of Anderson's paper is devoted to a derivation of (2.24). This 
same multiple integral, written in the form 

n+1 



(2.25) 



det 

l<'t,j<n 



ai\ 



1=1 



was evaluated at around the same time by A. Varchenko [152, 153] in his work on 
hyperplane arrangements. That (2.251 is equal to the integral in (2.241 is a simple 
consequence of the Vandermonde determinant, a fact made explicit in [129]. 
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Remarkably, in a 1998 paper by M.C. Bergere [20] proving a conjecture from the 



theory of Calogero-Sutherland models (see page 18 1 reference is made to (2.241 



citing neither Anderson nor Varchenko, but a paper of A.L. Dixon [37] published 



in 1905! Indeed, consulting [37] reveals both (2.24) — obtained via essentially the 



same analysis as that used in [3] — and its equivalent determinant form (2.25) 



A study by Forrester and E.M. Rains [57] provides additional links between the 
Selberg and the Dixon-Anderson integrals. These apply at the level of the corre- 
sponding normalized integrands, referred to as the Selberg and Dixon-Anderson 
densities. The former will be denoted by S„{a, /3, 7; t) for t = (ti, . . . , i„). 

The first point of note is that the computations of Dixon and Anderson can be 
interpreted as giving the density of zeros of the random rational function 

n+l 

Rn+i{x) := 



1=1 



a,; 



where the Wi are distributed according to the Dirichlct distribution — to be denoted 

^n+l[si, • • • , Sn+l] 

r(si + ■ ■ ■ + Sn+l) ''yv s,-l 

r(si)...r(5„+i) l\ ^ ' 

with wi, . . . , Wn+1 > such that wi + ■ ■ • + Wn+i = 1- Motivated by this interpre- 
tation, a family of random polynomials Aj(x), I < j < n were defined in [57] such 
that the zeros of Aj(x) have PDF Sj{aj, (3j,^;t) with aj := (n — j)^ + a, (3j :— 

{n — + p. Setting A^i{x) := 0, Ao{x) := 1, and specifying that {wq\wi\w2^) 
be distributed according to the Dirichlet distribution D^lPj, (j — 1)7, aj], the poly- 
nomials Aj{x) are determined by the random three-term recurrence 

(2.26) Aj{x) = wi^\x - l)Aj^i{x) + w^^^^'' xAj^i{x) + w[^^x{x - l)Aj^2{x). 

Let a t-^ ja/2 + 1, (3 7/3/2 + 1, and let the integrand of the Selberg integral 
be written in the form exp(^2^U) so that 

n o ^ 

u = -'^J2^ogt,-'^Y.^og\i-u\- iogi^^-*ji- 

Then in the limit 7 — > 00 the Selberg density crystallizes to the minimum of U 
subject to the constraint that < < 1 for each ti. According to a classical result 
of Stieltjes (see e.g., [6]) this minimum is unique and occurs at the zeros of the Jacobi 
polynomial Pn" ^'^ ^\x)- Indeed in the same limit the three term recurrence 



(2.26 1 is no longer random, and has solution ^j(x) = Pj" •''^^ (a;) with 



p(a,&)^2;) the Jacobi polynomial normalized to be monic [57]. 



The change of variables and limiting procedure giving rise to (1.8 1 reduces the 



Selberg density to the PDF (1.4). The Dixon- Anderson density permits a similar 



limit, and applied with ni-^n+l, = and ai — 1 this results in the PDF on 
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{t,}]+^ given by 



(2.27) 



1 



supported on 



/2^r(.si) • • • r(s„) ni<.<,<„(«. - %)^'+^^ 

n+1 n 

X 

^ j 



n-\-l n y ^ n-\-l i^ \ 

1=1 7 = 1 ^ 1 = 1 1 = 1 ^ 



ti > ai > t2 > ■ ■ ■ > an > tn+l- 

The corresponding limit of Rn+i{x) gives the random rational function 



(2.28) 



Rn+i{x) := X ~ Ha + y^^ 



where /ig has distribution N[0, 1] and fii the Gamma distribution T[si, 1]. Indeed, 



the fact that the zeros of (2.281 have PDF (2.271 can readily be checked directly 



by adopting the strategy of Dixon and Anderson, see [49,55]. Finally, the limiting 
form of the three term recurrence (2.26) is 



(2.29) 



C,{x) ^{x~ r)Cj-i{x) - s<^J-^^Cj-2{x) 



with C_i(x) — 0, Co(x) — I, r having distribution N[0, 1] and s^-'^ distribution 
r[j7, 1]. The random polynomial Cj{x) has as the PDF for its zeros the density 
(1.4 1 with n = j and /3 = 2j. 



It should be remarked that since (2.27) integrated over ti, . . . , tn+i gives unity, a 



limiting form of the Dixon- Anderson integral follows. Evans [49] used this, together 
with the strategy of Anderson, to give the first proof of the Mehta integral evaluation 



(1.6) which is independent of the Selberg integral. 



The random polynomial C„+i(x) can be interpreted as the characteristic poly- 
nomial for a family of random matrices defined inductively by [57] 



(2.30) 



n+l 



(») 



A. 



bi 



bn 



Here the A^"''' are the eigenvalues of M„, c has distribution N[0,1] and has dis 
tribution r[/3/2,l]. Indeed it is straightforward to show that the eigenvalues of 

bi and ai 



x\"-\ In the 

case (3—1, the invariance of members of the GOE with respect to conjugation by 



(2.30) are given by the zeros of (2.28), with /iq = c, /i 



orthogonal matrices shows that (2.30) is similar to GOE matrices, and an analo- 



gous understanding of the relationship between (2.30) in the case [3 — 2 and GUE 



matrices can be given. Moreover, it is generally true that a three-term recurrence 



Pj{x) = (a; - aj)pj_i{x) - bj_iPj_2{x) 
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with p-i{x) := 0, po{x) — I is satisfied by the characteristic polynomial for the 
tridiagonal matrix 

i-l fln-l bn-2 



&2 0,2 bi 

bi ai 

Hence Cn{x) is also the characteristic polynomial for the above random tridiagonal 
matrix with each Uj having distribution N[0,1] and with &| distributed as in ( |2.30 1. 
This is a result due to I. Dumitriu and A. Edelman [39], obtained without the use 
of (2.27). In this regard it should be mentioned that R. Killip and I. Nenciu [92], 
in a study which does not make use of the Dixon-Anderson integral (2.24), give 



the explicit construction of a family of random orthogonal matrices with eigenvalue 
PDF equal to the BC„ Selberg density, which itself is proportional to the integrand 
in (1.23). The methods of [39] and [92], which at a technical level proceed via 



a change of variables from a general tridiagonal matrix and unitary Hessenberg 
matrix to their eigenvalue/eigenvector decomposition, yield too the evaluations of 
the Mehta and Selberg integrals respectively. 



Dotsenko— Fateev integrals. In the course of studies in conformal field theory, 
V.S. Dotsenko and V.A. Fateev [38] were lead to consider the multiple integral 

p ^ ^ ^ n m n m 

J J J J Etn^-t^)'E-t\n-ifUIl(-^-t^r' 

[0,1]P [l,oo)"-P [0,1]'- [l.oo)™-'- i=lj = l 

X Y[ \U-tj\^^ Y[ |t, -r,fT'dti---di„dri---dr™, 

l<i<j<n l<i<j<m 

where PV denotes the principal value, a/a' — P/(3' — —7, 77' — 1 and < p < n, 
< r < m. Note that the case p = n and m = is, up to a shift by 1 in a and /3, 
precisely the Selberg integral. Dotsenko and Fateev evaluated the above integral as 
a product of gamma and sine functions reclaiming the Selberg integral as a special 
case. 

The method employed by Dotsenko and Fateev for evaluating their integral sug- 
gests an approach [54] to the Selberg integral by studying the simpler m — case 



(2.31) 5„,p(a,A7) := 




for < p < n. Note that S'„_„(q!, /?, 7) — S'„(q!, /3, 7), which is the Selberg integral. 
Also note that the change of variables ti 1-^ 1/ti for all 1 < i < ?i implies the 
transformation 

(2.32) 5„,p(a, P, 7) = Sn,n-p{l ~a~P-2{n- 1)7, /?, 7). 

Singling out the integration variable tp, viewing the integrand as an analytic 
function and replacing the interval [0, 1] by a contour along a positively oriented, 
indented semi-circle of infinite radius (with indentations at the branch points tp — 
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0, 1, ti, . . . , <p_i, tp+i, . . . , tn) yields the recurrence 

(2.33) Sn,p(a, P, 7) ^ — -. — — — 

n — p+l sin7rp7 sin7r(a + (p — l)7j 

X S'„,p_i(a,/3,7). 

Solving Sn,n in terms of Sn,o f^nd then using the transformation ( 2.32[ ) to eliminate 
Snfi in favour of Sn,n gives the following functional equation for the Selberg integral 

or o \ an a oi 1^ . -rr sin 7r(a + /? + (n + j - 1)7) 
S'„ a,/3,7 = Sn{l - a - /3- 2 71 - 1 7,/?, 7 II -. — — . 

The significance of this result is that it permits S'„, viewed as a function of a, to 
be analytically continued to a € C with the exclusion of the zeros of sin 7r(a + J7). 
Indeed, for sufficiently small but positive values of Re(/3) and {n — l)Re(7) the 
Selberg integral requires Re(Q;) > but the transformed integral permits Re(Q;) < 
1 — Re(/3) — 2{n — l)Re(7) which is greater than 0. As a consequence one can verify 
that 5„(a, /3, 7) when divided by the right-hand side of ( 1.1 1 is a bounded analytic 
function in the complex a-plane, and is thus independent of a. Symmetry then 
gives that this ratio is independent of /3 as well. That the dependence on 7 and n 
is correct is then verified by using 

lim a5„(Q;,/3,7) = 7iS'„_i(27 - l,/3,7), 

a fact already noted and used for the same purpose in Selberg's original proof [134]. 

In their paper Dotsenko and Fateev considered a further generalization of ( |1.1[ ) , 
referred to as the complex Selberg integral. This integral, which was also studied 
independently by K. Aomoto [8], can be written as an n-dimensional real integral 
with integration variables given by 2-dimensional vectors 

A„(a,/3,7) / TTlrlp^"-!) |u - rlpC^-D TT |r* - f.f drl • • • dr;, 

where u is an arbitrary unity vector. Dotsenko and Fateev as well as Aomoto showed 
that up to a product of trigonometric functions the complex Selberg integral factors 
as a product of two Selberg integrals 

2 / n s 1 -TT sin 7r(a + j^i) sin 7r(/3 + ^7) sin 7r( j + 1)7 

A„ a,/3,7 = S'„ a,^,7 — I I — , — — -— — , 

n! sin7r(Q; + p + (n + — 1)7) sin7r7 



provided (1.2 1 is supplemented by 

Re(a + /3 + (n - 1)7) < 1 and Re(a + /3 + 2(n - 1)7) < 1. 

K. Mimachi and M. Yoshida [111] (see also [110]) apply the theory of intersection 
numbers of twisted cycles to the conformal field theory study of Dotsenko and Fa- 
teev to give the evaluation of the product S'„(a, (3, 7)S'„(— a, —(3, —7), appropriately 
analytically continued. This is achieved without requiring the actual evaluation of 
the Selberg integral itself. 
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Jack polynomial theory. It has been known since the early 1970s [25, 148, 149] 
that (1.9 1 with (3 — 2^ — to be denoted exp(— 27!^) in analogy with (1.7) — is 
the absolute value squared of the ground-state wave function for the Schrodinger 
operator 

.-I — — > 

06^ 



" a2 1 

^ = -E7^+o7(7-l) 



E 



{0, - 6,)) ■ 



This operator, known as the Calogero-Sutherland Hamiltonian, describes a system 
of n identical quantum particles on the unit circle, with 6i G [0, 27r) for 1 < i < n 
the (angular) positions of the particles. The interaction between the particles is 
described by a 1/r^ two-body potential, 2|sin((6'i — dj)/2)\ being the cord-length 
between particles located at 9i and 6j. 

B. Sutherland [148] showed that the eigenvalue Eq corresponding to the ground- 
state wave function is given by Eq = n{n'^ — 1)7^/12. Subsequently he showed [149] 
that the conjugated operator 



(2.34) 



e^^{H-Eo)e 



-27 



E 



d 



where Xj :~ exp(i0j), admits a complete set of symmetric polynomial eigenfunctions 



P^^^'^\x). These polynomials, now referred to as Jack polynomials, depend on 



X — (xi, . . . , Xn) and are indexed by partitions A of at most n parts; A — (Ai, . . . , A„) 
with Ai > A2 > • • • > A„ > 0. With m\ denoting the monomial symmetric function 
indexed by A and < the dominance ordering on partitions, the Jack polynomials 
have the structure 



(2.35) 



P 



(1/7) 



{x) = mx{x) 



for some coefficients ax^ = a\^{'-^). 

One fundamental property of the Jack polynomials is that they are orthogonal 
with respect to the inner product 



(2.36) 



1 



(27r)' 



') n 

l<i<j <n 



-e'^^ |2^ d^i-.-d^^ 



where /(e'^) = /(e'^ 
norm evaluation let 



with {b)n = b{b+l) ■ 
be given by 



(6 



(2.37a) 
(2.37b) 



/(e"')ff(e- 

e'''"). To state the orthogonality as well as the quadratic 

\b]^x^ = \{{b+{l-^)l)x. 
h ri — 1) a Pochhammer symbol. Also let ca(7) and 0^(7) 

cx{l) = X{{a{s) 

CA(7)=n(«(^) 



^(5)7+7), 



sex 



where a(s) and l{s) are the arm-length and leg-length of the square s in the diagram 
of the partition A, and |A| is the total number of boxes in the diagram of A [102]. 
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Then 
(2.38) 



^p(l/7)^p(l/.)X ^^^^ 



c'a(7) 



r(l + n7) ^(1/^) 



[l + {n-l)j][^^ r»(l + 7) 



'(1"), 



where iJa^ is the Kronecker delta function and (1") is shorthand for (1,1,..., 1). 
The orthogonahty relation is consistent with, but not an immediate consequence of 
the operator (2.34) being self-adjoint with respect to (2.36). The complication is 
that not all eigenvalues of (2.34) are distinct. This degeneracy can be removed by 
introducing the mutually commuting Cherednik operators for I < i <n [27,40] 



Xt d 
7 dxi 



E 



(1 



(1 - Szi), 



j=i+l 



where Sij acts by permutation Xi and Xj and 1 represents the identity operator. 
Any symmetric combination of the ^i, and in particular nr=i(l ~ '^^^Oi the Jack 
polynomials as simultaneous eigenfunctions. 

The Cherednik operators can be used to construct the Jack polynomial shift op- 
erator — a special case of the shift operators studied by Heckman and Opdam, and 
used by the latter to prove the Macdonald integral and constant term conjectures. 
Properties of the Jack shift operator not only imply (1.12) or, equivalently, (1.13), 
but also the more general quadratic norm evaluation of the Jack polynomials cor- 
responding to ( |2.38| wi th A = ^ [84]. (For A = this yields \1.12\ .) With A{x) the 



7 



n(n-l)/2 



Vandermonde product (1.21) and Y± 
shift operators are defined by G+ :— A^^l + , G 

10 the 

5)7+1 



Y\i<i<j<ni^i - T 1), the Jack 
= y_A. They have an adjoint 



type property with respect to the inner product (2.36), 

(2.39) (G+/,5)7+i = (/,G_ff). 
Also, with 

(2.40) a±(7)= II (A,-A,±l-f(j-zTl)7) 



n 

l<i<j<7i 



and 6 the staircase partition (n — 1, 71 — 2, . 
Jack polynomials as 



(2.41a) 
(2.41b) 



G+P. 



(1/7) 



G_P. 



(1/(7+1)) 
A 



, 1, 0), the shift operators act on the 



(1/(7+1)) 



«A(7 + i)Pj; 
«A (7 + 1)^1+^7^ 



It follows from (2.39) and (2.41) that 



^pa/(7+i)) p(i/(7+i)) 



^ "a (7 + I) 
^+1 a+(7 + l) 



(1/7) 
\+s ■ 



P\ 



(l/7)\ 



A+i5 



and thus 



l/(7+fc)) D(l/(7+fc)) 



/ pVJ-/ \ rTA.;; pi, 
\^A ' 



\ _ /p(l/7) p(l/7)\ 

/7+fe ^ V\+kS ' ^A+fc<5 /^ 



fc-1 



n 



1 "a+j5 



(7 + ^ - j) 



Taking 7 = 0, using that p^""-* = mx (the monomial symmetric function) and 
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which is n! for fj, — X + k6, it follows that for nonnegative integer k 



(2.42) 



(i/k) p(i/fc)\ _ , 

A '-^A /fc - 



X+]k 



Using the evaluation formula [141] 



(2.43) 



P. 



(1/7) 



(1") 



j=0 "A+jfe 



ca(7) 



(fc-j) 



and the definitions (2.371 and ( 2.40[ ) it is now a straightforward exercise to verify 
that for 7 = /c (2.38) coincides with (2.421. Analytic continuation off the integers 
is then required to establish (2.381 for all Re(7) > — 1/n. 

A further fundamental property of the Jack polynomials is R.P. Stanley's Cauchy 
identity [141] 



/ \ n m 



(2.44) 

-• — ■ C, 1 'V I 

-i=lj=l 

where x = {xi, . . . ,x„), y — (yi, . . . ,ym,)- Setting y — (1™), using the evaluation 
formula (2.431 with n i— *■ m and a standard analytic argument to replace m by aa, 
leads to Z. Yan's binomial theorem for Jack polynomials [161] 

r i(") " 1 

(2-45) V ^ Pt\x) = n j^^. 

^ ' ^ c\(a) ^ ^ ' (1-Xi)" 



This, together with the orthogonality (2.381, the property 

^S+a,...,A„+a)(^) = (^1 • • • ^nrPt\^h a = 0, 1, 2, . . . 

and the gamma reflection formula, implies a generalization of the Morris integral 
dTlTt [17] 

1 



(2^) 



f ... f pf/^)(_e'^) J]e5'«-(-'')|l+e'^T+' 



1 = 1 

X n |e'^'-e'^^f^d0i---d0„ 

l<i<j<n 



[l + a+(n-l)7]l:^ 



Pf/^^(nM„(a,&,7). 



Applying (1.15) finally results in a generalization of the Selberg integral 

(2.46) ■■■ pi'^'\t)i[tr\i-kr-' n i^.~^.rdii-..dt„ 

•^0 -^0 »=1 l<i<J<n 

= ^ ["+("7^^^Jf',, pf/-Hn5„(a,A7). 

[a + /3 + 2(n-l)7]i^) 

This evaluation is usually referred to as Kadell's integral [83] after its first prover, 
but as a conjecture is due to Macdonald [101, Conjecture (C5)]. When A = (l*"), 
in which case the Jack polynomial is nothing but the rth elementary symmetric 
function, the above is known as Aomoto's integral, who used it to give what is 
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arguably the first elementary proof of the Selberg integral [7]. A proof of Kadell's 
integral along the lines of Anderson's proof of the Selberg integral as described on 
page 13 has recently been obtained in [156] through use of the Okounkov-Olshanski 
integral formula for Jack polynomials [116] 

l\ r(A. + (n-*)7)r(7) ^JJ-^^ 

l<i<j<n-l 2=1 j = l 

where Y denotes the domain xi < ?/i < a;2 < ■ • • < < j/n-i < Xn and A is a 

partition of at most n — 1 parts. 

An open problem, settled only in the 2- variable case [50], is to find (and prove) 
a finite field analogue of Kadell's integral. 



The binomial theorem for Jack polynomials (2.451 can succinctly be written in 
hypergeometric notation as 



i-Fo 



(7) 



n 



where i-Fq^^ is an example of a hypergeometric function with Jack polynomial 
argument 

Ui,...,^/"J-^[6,]i'')...[a'') c',(7) ■ 

When n = 1, so that a; is a scalar, this function reduces to the classical hypergeo- 
metric function rFg- For general n, hypergeometric functions of the type (2.47) have 
their genesis in the work of A.G. Constantine [30], C.S. Herz [72] and R.J. Muir- 
head [114], but were first studied in their full form presented here by Kaneko [85], 
A. Koranyi [94] and Yan [161]. The case r = 2, s = 1 of (2.471 shares many prop- 
erties in common with its n = 1 counterpart, the Gauss hypergeometric function. 

One such property is that 2-Fi*'^'' solves the n-dimensional analogue of Euler's 
hypergeometric equation. Specifically, Yan [161] and Kaneko [85] independently 
showed that 2Fi-^\a,b;c;x) is the unique symmetric function, analytic at the ori- 
gin, that solves the system of n partial differential equations 

d'^F ( \ dF 

Xi{l - + (^c - (n - 1)7 - (a + 5 + 1 - (n - l)'^)x.ij — - abF 

V- 1 ^dF ^dF 

for 1 < i < 71. 

In the one-variable theory the Gauss hypergeometric function admits an integral 
representation due to Euler [45] 



22 



PETER J. FORRESTER AND S. OLE WARNAAR 



for Rc(c) > Rc(6) > with a branch cut in the complex z-plane from 1 to infinity. 
When z = 1 the integral on the right is the beta integral ( 1.3 1, resulting in a closed 



form evaluation of the 2F1 due to Gauss. In the multivariable theory an analogous 
results holds, where now the key integral-evaluation is provided by the Selberg 
integral. Multiplying both sides of Kadell's integral by z^^^[a]'"^^ /c'^{-y), summing 



the left-hand side using the binomial theorem (2.45 1, and using the definition (2.47) 



on the right-hand side, shows that Euler's integral extends to 



(2.49) 2Fi(-)(":';(z"))^n 



r(c-j7)r(i + 7) 



T{b - n)T{c - b - j7)r(i + (.7 + 1)7) 



„1 „1 n 

Jo Jo .I 1^.^.^., 



l<i<j<7l 



with a = b — {n — 1)7 and (3 ~ c — b — (n — 1)7. Evaluating the z = 1 instance of 
the integral by the Selberg integral (which, incidentally, follows by taking z = in 
(2.49 1) implies a generalized Gauss summation [161] 



^p^h) ^. (i")"! = Yl ^ ■^'^)^(^ -a-b-jj) 



Ln r(c - a - j7)r(c - 6 - j7) 



In addition to relating to the Selberg integral, the Cauchy identity also gives 
rise to a special limiting case of the Selberg density, referred to as the Laguerre 
PDF [57,77]. To motivate the origin of this we remark that the Jack polynomial at 
7 = 1 is equal to the Schur polynomial s\ while ca(1)/c';^(1) = 1. The normalised 
summand of (2.441 then reads 



sx{x)sx{y) YlYlC^ - X.Vj) 
i=ij=i 

which may be recognised as the measure on partitions induced by the Robinson- 
Schensted-Knuth correspondence, see e.g., [78,118]. As such the Schur measure 
holds a special place in certain studies relating to the representation theory of the 
symmetric group [24, 127]. 



For general 7 the normalised summand of (2.441 implies the more general mea- 
sure on partitions 



(2.50) 



ca(7) 



P 



(1/7) 



(1/7) 



(2/) n 11(1 



where n < m and A a partition of at most n parts. To obtain the Laguerre PDF 
[57,77] one needs to specialize x and y in (2.501 to 



1 < i < n and yj — <?^^^, 



I < j < m. 



Use of the Jack polynomial evaluation formula ( 2.43 1 allows all terms in ( 2.50 ) to be 



made explicit. The remaining step is to take the scaling limit, turning the discrete 
measure on partitions into a continuous one on the positive real line. This is done 
by setting q = exp(— 1/L), introducing the scaled variables tj := Xj/L and then 
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taking the large L limit for fixed tj. One finds that (2.501 multiplied by L" tends 
to the PDF, supported on ti > i2 > • • • > t„ > 0, 

n 

T4^„((m-n+l)7,7)i^ i<,<.<„ 



where Wn is the normalization 
(2.52) Wn{a,^)=-\{ 



" 'r(a + j7)r((j + l)7) 



n\ J-J- r(7) 

To obtain this same result as a limit of the Selberg density, order the integration 
variables in the latter and write /3 — L. Then the change of variables ti i— > ti/L 



followed by the limit L —>■ oo results in ( 2.51 ) after identifying a with {m — n + 1)7 



This limiting "Laguerre" case of the Selberg integral, leading to the evaluation 
( 2.52| , is contained as equation (1) in the letter from Selberg to Dyson reprinted 



on page[5j 

Askey and D. Richards [13] (see also [107]) have shown that after some fairly 
straightforward manipulations and a change of variables the Laguerre limit of the 
Selberg integral leads to the evaluation 

(2.53) / l[tr'{l-Y.^^) n \U-t,\^''dh---dt^ 



i—1 ^^'i<j<n 



r(/3) ^r(a + .77)r(l + (j + l)7) 



n 



T{an + f3 + n{n-l)-f) ^J-^ r(l + 7) 

where D is the domain ti > for 1 < i < n such that ti + ■ ■ ■ + tn < 1, and 
Re(a) > 0, Re(/3) > 0, Re(7) > - min{l/n, Re(a)/(n - 1)}. 



According to Askey and Richards, the first statement of (2.531 is due to Selberg at 



a meeting held in Sri Lanka during December 1987. The intriguing point is that 



while Selberg did not give his derivation of (2.531, he is reported to have said that 
it was different to ( |1.1[ ), and had the advantage of working in the finite field case. 
The derivation given in [13] does not work in the finite field case, and therefore 
must be different to that known to Selberg. 

q-Integrals and constant terms. Motivated by the Selberg integral and its suc- 
cess in dealing with Dyson and Macdonald type constant term identities, Askey in 
1980, was led to consider several (/-analogues of the Selberg integral and to study 
connections to g-constant term identities. In fact, one learns from [12] that he had 
earlier spent time searching for a proof of the Mehta integral upon its appearance in 
the problem section of SI AM Review. To describe some of Askey 's work we require 
the multiple Jackson or g-integral 



„i „i 00 
/ •••/ /(i)d,ti...d,i„:=(l-g)" fil')l''' 

"'O Jo ,., .. _n 



iH hfc„ 

fei,...,fc„=0 

with t = (ti, . . . ,t„), q'^ = (q'^^, . . . jq'^") and < <; < 1, and where it is assumed 
that the multiple sum on the right is absolutely convergent. Also needed is the 
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Q'-shifted factorial 



{a;q)z = n Y 



1 — aq^ 



3=0 



aq 



for z E C Probably the most important of the g-Selberg integrals considered by 
Askey is 

(2.54) 5„(a,/3,7;g) 

Jo Jo l<i<j<n 

It is immediate, at least formally, that 

lim S'„(a,/3,7; g) = 5'„(q;,/3,7). 

For Re (a) > 0, 7 a nonnegative integer, say k, and /3 G C excluding the nonpositive 
integers Askey conjectured (and proved for k — 2) that [11, Conjecture 1] 



S.ia, A k; q) = ,-(S)+-^(S) n + 0- - mVqW + + ^-fc) , 

j=o r,(a + /3+ (n + J - 2)fc)rq(l + fc) 

where Tq{x) is the g-gamma function 

_ (g;g).-i 
''^ ^" (l-g)-i- 

For Askey's other g-Selberg integrals and many further results relating to Jackson- 
integral type extensions of beta integrals, see [5, 9-11, 48, 76, 79, 86, 155]. 

In 1988 Askey's conjecture was proved independently by L. Habsieger [70] and 
K. Kadell [80]. Both then applied the q-analogue of the identity (1.15) to (2.54) to 
obtain a generalization of the Morris integral (1.17). Expressing this integral as 
a constant term identity they thus proved Morris' g-constant term conjecture [113] 



W{xf,q)a{q/xf,q)b ]J {x^/xj■,q)k{qxJ/xi■,q) 

n-1 -p 

r„(i 



T,{l + a + b + jk)Tg{l + {j + l)k) 
^.^0 . ,(1 + a + jk)Tg{l + b + jk)Tgil + k 



When a = & = this is precisely the A„_i case of the g-Macdonald constant term 
conjecture [100, Conjecture 3.1] 



(2.55a) 
where 



fc I 
CT n ll{l-q'-'e-''){l~q'en^l[ 



dik 
k 



l-q 



is a g-binomial coefficient, and <i> is a reduced (finite) root system. To also include 
the root systems of type BC one again needs the numbers s{a) as defined above 
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equation (|1.22[) [100, Conjecture 3.4] 



(2.55b) 



CT n 11(1 

a£^+ 1=1 



I r 



«)(1 



(■i-l)s(a) + l a\ _ 



dik 
k 



The A„_i case of (2.55), was in fact proved prior to the work of Habsieger 
and Kadell by Zeilberger and D.M. Bressoud [162], who proved the more general 
q-Dyson conjecture 



CT n 

l<i<j<n 



(9;9)ai •••(<?; g)a„ 



formulated by G.E. Andrews [4]. 

R.A. Gustafson [67], at around the same time as Anderson's work on the Dixon- 
Anderson integral, made use of a further ^-generalization of (2.24) and invented the 
same general strategy as used in [3] to derive the BC„-type constant term identity 



U) = 2"n! W 



\ {t^\q)ao{q\q)ooY{i<r<s<ii^^ '^trts]q)c 



1 Or— 1 (^''^'^^ ' ^)oo(^r*^2 ' 9)00 



(2.56) CTA{x;t,ti, 
where 

n 

(2.57) A{x;t,tu.--,h)-^l[ 

l<i<j<n ^^'^^'^j i9)oo(tej Xj , q^oaitXiXj , q)oo{tXi Xj , q) oo 

This is a generalization of the so-called Macdonald-Morris constant term identity 
and implies the B„, C„, D„ and BC„ cases of (2.55) through specialisation [67]. 

Most other cases of (2.55) were proved on a case by case basis, often using 
methods based on g-integrals of Selberg type [48,61,62,69,82,145,163-165], but the 
three exceptional root systems Eg, E7 and Eg refused to surrender until Cherednik 
gave a uniform proof for all reduced root systems based on his theory of double 
affine Hecke algebras [28]. 

Returning to Askey's g-Selberg integral we remark that the density function 
corresponding to the integrand of (2.54) can be deduced from Macdonald polyno- 
mial theory [102], following a procedure similar to that of deducing (2.51) from 
(2.501 [57]. Macdonald polynomials P\{x;q,t) are generalizations of Jack polyno- 
mials — the latter being reclaimed according to lim^^^i P\{x; q, q'^) — pI^^^''\x) — 

■d9i- ■ -dOn cx S\^. 



exhibiting the structure (2.35) and the orthogonality 

(e-(^'-^^);g)o< 
(te^(»'-«.);g)c 



Px{e^';q,t)P^{e'^';q,t) J] 



l<i<j<.n 

It is the connection between afhne Hecke algebras and Macdonald type orthogo- 
nal polynomials that is at the heart of Cherednik's proof of the g-constant term 
conjectures for arbitrary root systems [28], see also [103]. 
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Multi variable orthogonal polynomials. By the change of variables t ~ {1 
x)/2 and a shift in a and /? by 1 the Euler beta integral (1.3) takes the form 



J(a,/3) 



The integrand on the left is the weight function of the Jacobi polynomials Pn°''^\x) 
[6,75]. Up to normalization these are the unique functions, analytic around the 
origin, solving the second order ODE 



(1 - x^)y" + (^l3-a- x{a + P + 2)jy' + n{n + a + l3 + l)y 
Defining the inner product 



= 0. 



(/,.9)a,/3 := J J{x)g{x){l-xni + xfdx 
the Jacobi polynomials satisfy the orthogonality relation 

All other classical orthogonal polynomials, such as the Laguerre and Hermite poly- 
nomials (corresponding to weights a;"exp(— x) and exp(— 2;^) respectively) follow 
from the Jacobi polynomials by taking appropriate limits. 

Several people have studied multivariable generalizations of the Jacobi, Laguerre 
and Hermite polynomials [16,33,35,41,95-97, 104, 121]. The most general of these 
are the multivariable Jacobi polynomials P^^^'^'^^^x) which arise as the eigenfunc- 
tions of the operator 



E((l-D^ + (/?-«-..(a + /? + 2)) )+2,5: 



1-xl d 



The Pi^^^^^^x) are orthogonal with respect to an inner product with weight func- 
tion derived from the Selberg integral. With 
(2.58) 

/n 
f{x)9{x)'[[{l-X^)°'{l+X^)'^ Y[ \X^- Xjl'^'' dxi- ■ - dXn 

[_X l]n l<i<j<n 

for X = (xi, . . . , Xn), the multivariable Jacobi polynomials satisfy 

The quadratic norm evaluation can be computed explicitly in term of Pochhammer 
symbols and gamma functions using the shift operators of Heckman and Opdam 
[119]. From the Selberg integral it of course immediately follows that 

(1, 1)„,^,^ = 2"("+''+i+("-i)^)5„(a -I- 1, ;3 -I- 1, 7). 



Two important limiting cases of the inner product (2.581 are 

/n 
/(x)5r(a;) J|e-^- H l^i ~ Xj\^'' dxi ■ ■ ■ dxn 



l<i<j<7i 
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and 



/n 



rn \7i l<i<j<n 

The corresponding families of orthogonal polynomials are referred to as the multi- 
variable Hermite and Laguerre polynomials, respectively. In particular note that 



(1,1)^ = 2(3)^W2f„(7) 



with Fn the Mehta integral (1.5 1, and 



= n! VF„(a+ 1,7), 



with Wn given by (2.521 



All of the orthogonal polynomials mentioned above admit g-analogues. In the q- 



theory the role of the Selberg integral is played by Askey's g-Selberg integral (2.54) 
and generalisations thereof. In the case of the Jacobi polynomials these g-analogues 
are known as the multivariable big and little g-Jacobi polynomials, and were intro- 
duced by J.V. Stokman [146]. Stokman [147] also showed how the big and little 
g- Jacobi polynomials arise as special limits of the Koornwinder polynomials [93]. 
The latter are multivariable analogues of the Askey- Wilson orthogonal polynomi- 
als [14] and may be viewed as the generalizations of the Macdonald polynomials to 
the root system BC„. The relevant inner product in this case is given by [34,93] 
(2.59) 

{f.9)tM,-M-= 77^ I ■■■ I /(e'^)5(e'")A(e'-;t,ti,...,i4)dxi---dx, 



(27r) 

where exp(ia::) = (exp(ia;i), . . . , exp(ia;„)), A(x; t,ti, . . . ,14) is the weight function 



(2.571 of Gustafson's constant term identity, and / and g are BC„ symmetric func- 
tions. {f{z) is BC„ symmetric if /(exp(ix)) is symmetric under signed permuta- 
tions of a: = (xi, . . . ,Xn)-) The evaluation of (1, \)t,....ti is of course provided by 



the right-hand side of (2.56) 



3. Recent and current research directions 

The case of 7 a positive integer. Two recent studies have identified special 
features of the Selberg integral for 7 a positive integer. The first of these is due 
to J.-G. Luque and J.-Y. Thibon [98], and exhibits an inter-relation with a special 
class of hyper determinants. The second, due to Stanley [143], gives a probabilistic 
interpretation of Selberg 's integral. 

For a fcth order tensor A = [yliji2...ij.] on an n-dimensional space (so that 1 < 
ip < n) the hyperdeterminant has been defined by Cay ley (see references in [98]) as 

n 

detfc(A) := ^ £(0-2) • • • e(CTfc) J| A,^<,^(,)^...^^^(i) 

where e{a) denotes the signature of the permutation a. For k odd this vanishes 
while k — 2 corresponds to the usual definition of a determinant. 

For an arbitrary measure /i(x) it is easy to see by use of the Vandermonde 
determinant formula that in the so-called Hankel case 



^ili2---ik Mii+i2H hifc— fcj ■ / X'^ dfi(x) 
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the corresponding hyperdeterminant is equal to a multiple integral, 



det2fc(A) = 




For d/i(x) — x°' ^(1 — x)^ ^dx on a; G [0, 1] this is precisely the Selberg integral 
with 7 = fc an integer. 

The probabilistic interpretation of the Selberg integral given in [143] applies for 
a = /5 = 1 and 7 a positive integer. In fact, as communicated to us by Stanley 
[144], this same probabilistic interpretation, extended to a,/3 general nonnegative 
integers, is already implied by appropriately interpreting supplementary problem 
1.25 of [142]. Following [144], the setting is to choose labelled points independently 
and with uniform distribution from the interval [0,1]. Specifically, for each 1 < 
p < n and such that I < i < j < n, choose a — I points labelled j/p, (3—1 
points labelled Zp, n points labelled t and 27 points labelled Uij. Let ti be the i-th 
smallest point labelled t. The probability that any one of the points labelled Up 
is to the left of ti is 1 — ti; the probability that any one of the points labelled aij 
is between ti and tj for z < j is — tj. It follows immediately that the Selberg 
integral S'„(a,/3,7) is the probability, Ps say, that all the points labelled t/p are to 
the left of tp, and all the points labelled Zp are to the right of tp, and all the points 
labelled Oy- lie between ti and tj. Note that this statement remains valid for 27 
an odd integer. An equivalent formulation (the one given in [142]) is to regard the 
selection of the labelled points uniformly and independently from [0, 1] as a random 
re-arrangement of the symbols themselves (i.e. the t's, j/p's, ^p's and a^j's), in which 
case Ps corresponds to the probability that the re-arrangement complies with the 
prescribed rule. 



Random matrix theory. A number of interplays between random matrix theory 
and the Selberg integral appearing in papers published in the last few years were 
discussed previously under the heading of the Dixon-Anderson integral. Below two 
further applications of the Selberg integral to random matrix theory as they have 
occurred in current works will be outlined. 

The first of these is a study by Forrester and Rains [58] focusing on the family 
of multi-dimensional integrals 

(3.60) /„,p(x) 

•= / / f\tr\^"t'^f~'\^-t^r' IT \u-t,\'^dh---dt^ 



for < p < n and x G [0, 1]. In the case r = 1 the integral (3.60 1 is proportional to 
the probability that for the Selberg density, interpreted as an eigenvalue PDF, there 
are p eigenvalues in the interval [0, x] and n — p eigenvalues in the interval [x, 1]; for 



T = 1 -I- 27 the integral (3.60) relates to the derivative of this quantity. The case 
T = 27 = 1 of this was first studied in the mathematical statistics literature for its 
relevance to canonical correlation analysis [32] . 

Theory connecting /„_p to a certain Fuchsian differential equation [32,53,109] im- 
plies that the integral is expressible as a linear combination of Frobenius solutions. 
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These are solutions to the differential equation of the form 

oo 

gi{x) = x"' ^ Oi^fe x^, = i(a - 1 + T + (i - 1)7) 

for < * < and are normalized such that 
(3.61) 



lim = 1 for Re((7,) > 0. 

2^—0 In,p{x) 



A basic task — essentially equivalent to finding the monodromy matrix for the 



basis of integral solutions of the matrix Fuchsian system, of which (3.60) forms the 
top row — is to give the explicit form of the coefficients Cp^i in the expansion 



(3.62) 



^7i,p{x^ — ^ ^ Cp^i (/i(x). 



One approach to this problem is to seek a regime in parameter space such that 
for a; — > the leading behaviour of In,p{x) is proportional to x"^ . This is achieved 
by changing variables tj = xuj for 1 < j < i where p + 1 < i < n. A simple scaling 
of the integrand then shows that 

In,p{x) ^ x'^'Sn~i{a + T + 27 - 1 , /?, 7) 5*,;,^ (a, T, 7) , 



where Sn^i is the Selberg integral and 5";^^ the Dotsenko-Fateev integral (2.311 



Recalling the normalization (3.61) and the recurrence (2.33) allows the sought co- 
efficients to be calculated as 



0, 



l — p 



sin 7r(z — J + 1)7 sin ir^a + {i — 
jfj sin 7rj7 sin Tr(a + r — 1 + (2i — j — 1)7) ' 

thus solving the problem at hand. 



< i <p- 1 
p < i < n 



The second of the applications stems from a question posed by Balint Virag at a 
recent AMS-IMS-SIAM summer research conference [1] . By way of background to 
his question, let us recall a result of Mehta and Dyson [108] which gives the circular 
ensembles identity alt(C0E2„) = CSE„. Here alt is the operation of integrating out 
every second eigenvalue, and the subscripts on the names of the ensembles indicate 
the total number of eigenvalues. Let us more generally introduce the notation CE^.„ 
for the PDF ( 1.9 1. The question posed by Virag was to investigate extensions of the 



result of Mehta and Dyson, in which blocks of eigenvalues in CE^ „ are integrated 
out to obtain another circular ensemble CE^' „/. The Selberg integral is relevant 
for this purpose. 

Let p{k; s; CE^_„) denote the PDF for the spacing between eigenvalues which are 
(fc + l)-st neighbours in the ensemble CE^^^v- Let alt„i(CE^^„) denote the joint 
marginal distribution of every m-th eigenvalue in CE^^„. With this notation, if it 
were true that 

(3.63) altm(CE^_„„) — CE^'^„ 
for some m,/3, /?', then 

(3.64) p{mk + m - 1; s; CEp^mn) = Pik; s; CE^/_„). 
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Now the fc-point correlation pk is obtained from (1.9 1 by integrating out the 
variables Ok+i, . . . ,9n and multiplying by n!/(n — fc)!. It follows from this definition 
that pk is related to the small s expansion of p according to 

(3.65) pik; s; CE^,„) ^ ^ ' ' ' ^ Pfc+2(0, s, 0^, . . 

But for di, . . . ,9n small the definition of pk also implies that 

1 n\ M„,fc(fc/3/2,fc/3/2,/3/2) 
(2^)fc(n-fc)! M„(0,0,/3/2) 



(3.66) pk{ 







o)d(?i 



■d0k 



n I 

l<i<j<k 



1/3 



where M„ refers to the Morris integral (1.171. Substituting (3.661 in (3.651 and 



scaling the integrand an example of the Selberg integral is obtained, giving the 
formula 



(3.67) p(fc;s;CE^,„) 



k+l3{k+2){k+l)/2 



1 {n-iy. 

(27r)fe+i fc!(n-fc-2)! " 
M„_fc_2((fc + 2)/3/2, (fc + 2)/3/2, /3/2) 
M„(0,0,/?/2) 



Skip +1,13 +1,13/2). 



Using the gamma function evaluations (1.1 1 and (1.181, together with the dupli- 
cation formula for the gamma function, one can check that in the case m — r + I, 



13 = 2/{r + 1) and f3' = 2{r + 1), ( |3.64| is compatible with ( |3.67[ ). Thus, this in- 
vestigation based on the Selberg integral reveals parameters for which the validity 
of ( |3.63[ ) may be expected. One can in fact proceed further and prove, using a 
generalization of the Dixon- Anderson integral, that for these parameters (3.631 is 
indeed valid [56]. 

KZ equations and the Mukhin Varchenko conjecture. On page[22]we have 
seen that hypergeometric integrals of Selberg type arise naturally as solutions of 
(systems) of partial diff'erential equations. There is a well-developed theory ex- 
tending much of this to the setting of partial differential equations — referred to as 
Knizhnik-Zamolodchikov (KZ) equations — based on Lie algebras [43,115,132,154]. 

Let be a simple Lie algebra of rank n, with simple roots, fundamental weights 
and Chevalley generators given by a^, and ei,fi,hi for 1 < i < n. Let Vx and 
Vf^ be highest weight representations of g with highest weights A and p, and let 
u = u{z, w) be a function taking values in V\ (g) solving the KZ equation 



du 
dz 



w 



du 
dw 



n 



u, 



where S 0(8)0 is the Casimir element. (For the sake of simplicity we only consider 
the KZ equation in two variables, z and w; for the more general case of p variables 
zi,...,Zp, see e.g., [154].) Let Sing;^^^[:^] denote the space of singular vectors of 
weight in Vx(^ Vfj, 

Sing;^^^[i/] := {v eVx(g)V^ : h.^v = v{h.i)v, e^w = 0, l<i< n}. 

Then, according to a theorem of V.V. Schechtman and Varchenko [132], solutions 
u with values in Sing^;^ ^[X + p — ^jO^i] £^re expressible in terms of multiple 

hypergeometric integrals 

u(z, w) = ^ ULj{z, w) f^vx (g) /■■'v^ 
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with coordinate functions u/j given by 

uij{z,w)= / <^^/''{z,w;t)Aij{z,w;t)dti---dtk. 
Ji 

Here k := fci + • • • + fc„. t := (ti, . . . , tk), the sum is over all ordered multisets 
/ and J with elements taken from {1, . . . ,n} such that their union contains the 
number i exactly fcj times, vx and t;^ are the highest weight vectors of V\ and V^, 
f^v := (Iljg/ ,fi)v and 7 is a suitable integration domain. 

The functions $ and Ajj in the integrand of ujj are explicitly known. Ajj is 
a rational function whose general form is too involved to explicitly state here (an 
example will be given below), and the function $, known as the master function, is 
defined as follows. The first ki integration variables are attached to the simple root 
ai, the next ^2 integration variables are attached to the simple root a2, and so on, 
such that at^ '■= en if fci + h ki-\ < j < ki + ■ ■ ■ + kt. With this understood 

k 

<P{z,w; t) = {z- J[{ti - z)-^^'°''i\ti - 

i=l 

X n (^-^.•)'"'*'"*^^ 

l<i<j<k 

where ( , ) is the bilinear symmetric form on f)* (the dual of the Cartan subalgebra 
[)) normalised such that (6, 9) = 2 for the maximal root 9. 

The simplest possible example of a KZ solution contained in the Schechtman- 
Varchenko theorem corresponds to the rank 1 Lie algebra g = 5(2 = Ai, with 
simple root a\ and fundamental weight Ai = a\/2. Taking A = miAi and ^, = 
m2^i it follows that ti{z,w) takes values in the space of singular vectors of weight 
(mi + m2 - 2/ei)Ai. Since n = 1 it follows that / = {V} and J = {I'^i"''} with 
1 < r < fci, so that w/j, Au and can simply be denoted by Ur, Ar and /''. (In 
the case of rank one there is no need for the index in /i). Using n instead of ki (so 
that n no longer denotes the rank of the Lie algebra) and writing vi and V2 instead 
of v\ = VmiAi and = WmjAa o^e finds, upon the assumption that z <w are both 
real, 

n 

U{Z, w) = ^ Ur{z, W) fvi O /"~''W2 
r=0 

with 

/n 
Ar{z,w; t) ^{ti - z)-"'''''{ti - w)-'^^/" 

X n iti-tjf''"dU---dtn. 

l<i<j<n 

Here the domain of integration is the simplex 7 = {<; e M"| z < tn < ■ ■ ■ < ti < w} , 
and the rational function Ar{z,w; t) is given by 

Ar{z,w;t)= n^n^- 

|7|=r 

The coordinate functions are easily recognised as generalizations of the Selberg 
integral. In fact, for the extremal cases r = and r = n they are exactly the 



32 



PETER J. FORRESTER AND S. OLE WARNAAR 



Selberg integral. When r — 0, for instance, 

p n 



T i=l 

l<i<j<n 



Making the change of variables ti = {w ~ z)si + z for 1 < i < n this yields 

{-1)^{Z-W)^ / mi TO2 1 
Uo{z,W) = 1 Sn[l , , - 



where A = n(n — 1 — mi)/ k + n and B — (mim2 — 2n{mi + + 2n{n — 1))/(2k). 

In 2000 E. Mukhin and Varchenko [115] formulated a surprising conjecture re- 
garding the scaled master function 

They conjectured that if the space Sing[A + — J27=i ^i^i] of singular vectors is 
one-dimensional, then 

(3.68) [\m\'^^dh---dtk 



is expressible as a product of gamma functions. Neither the exact integration 
domain nor the specific form for the product of gamma functions is contained in 
the Mukhin- Varchenko conjecture. 

For 2 — 5I2 ~ Ai the conjecture corresponds to the evaluation of the Selberg 
integral. For g = sl„+i = A„, and Sing;^^^[A + - Z^Li "^i*^ X = Ai, fi = 
Y17=i the conjecture simply follows by iterating the beta integral (1.3 1, see 
[115]. For Q = B„, C„ or D„ and 

n-2 I (2, 2) B„ 

SingAi 2Ai - rQ;„„i - sa„ - ^ a^j with (r, s) ^ I (2, 1) C„ 

'=1 i(l,l) D„ 

(corresponding to the tensor product of the vector representation of g) Mimachi 
and T. Takamuki [112] established the Mukhin- Varchenko conjecture iterating the 
Selberg integral for n = 2 (B„ case) or the beta integral (C„ and D„ cases). 

In 2003 V. Tarasov and Varchenko employed KZ equations and the closely re- 
lated dynamical equations to settle the conjecture for g = 3(3 A2. In recent work 
by Warnaar [157,158] an approach to the s[„+i = A„ case of the Mukhin- Varchenko 
conjecture was developed, based on the theory of Macdonald polynomials and gen- 
eralized hypergeometric series. Specifically, the integral (3.681 for g = A„ can be 
evaluated in closed form when A = A„A„ and /i = /XjAi (or when A — AiAi and 
/i = Stripping the integral from its Lie algebra notation and using 

and (3i (1 < i < n) for exponents in the integral (so that the no longer denote 
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the simple roots) the A„ Selberg integral can be stated explicitly as 



n-l 



4=1 



S = l 



n n 



r(/3, + --- + f3r + ii + s-r- 1)7) 



l<s<r<n i 
n k 

X 



nn 



k 



r+l 



2)7) 



r(7) 



Here fci, . . . , fc„_|_i are nonnegative integers such that = and ki < k2 < ■ ■ ■ < 
kn, the exponents ai, q;„, /3i, ...,/?„, 7 G C such that ai = ■ ■ ■ = = 1 

and such that both sides of the identity are well-defined. Furthermore, i'*^ = 
{t^i \ . . . , i^*'') is the set of variables attached to the sth simple root of A„, 



A(u) 



{u^-Uj) and A(m, w) = ]^ - Wj) 

i<i<j<;„ i=ij=i 

for sets of variables u — (ui, . . . ,uij) and v = {vi, . . . ,vi^), and dt = d^^^^ • • • di^") 
with dt^'*^ = dt[^^ ■ ■ ■ dt^'j^^ so that the integral is (fci + • • • + fc„)-dimensional. 

Not yet specified in the A„ Selberg integral is the domain of integration, which, 
unfortunately, is rather involved. A key ingredient is the set of maps 

Af, :{l,...,fcj^{l,...,fc,+i} 

such that 

Ms{i) < Ms{i + 1) and 1 < Ms{i) < fcs+i - kg + i. 

A standard counting argument shows that there are exactly Ck^^-^^k^ admissible 
Ms, where Cn.k is the row (n,k) entry in the Catalan triangle, or, equivalently, 
the number of standard Young tableaux of shape (n, k). Given Ms fix an ordering 
among the t'f'' and tj''^^'^ as 



(3.69) 



where :~ 00. Given admissible maps Mi,...,M„_i define -Djvfi 

the set of points 



as 



(1) 



.(1) .(2) 



.(2) 



An) 



i1) 



, Ul , . . . , Ij,^ 



such that (3.69 1 holds for all 1 < s < 71 — 1 and 

< 4!? < • • • < ^ < 1 

holds for all 1 < s < n. Then the domain of integration, written as a chain, is given 

by 



-'Mi,...,M, 

Mi,...,M„_i 



iWj^Mi,...,M„_i' 



where 



n— 1 /i: 



Mi,...,M„. 



.(7)=nn 



^ sin(7r(« + /c. 



s+l 



Af.W + l)7) 



sin(7r(i + fc^+i - ks)-f) 
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In complete analogy with the ordinary Selberg integral, the evaluation of the A„ 
Selberg integral can be generalized to include a Jack polynomial in the integrand, 



thus generalizing the Kadell integral (2.461, see [158] 



Elliptic Selberg integrals. In the last few years there has been rapid progress 
in the field of elliptic generahzations of hypergeometric series, see [63, 140]. Classi- 
cal hypergeometric series X)^o characterized by the ratio c„+i/c„ being a 
rational function of n. Their elliptic counterparts have the same ratio equal to an 
elliptic function of n. 

As well as the classical hypergeometric series permitting elliptic generalizations. 



so do related integrals such as the Euler beta integral (1.3 1. In the elliptic theory 
the ordinary gamma function must be replaced by what is known as the elliptic 
gamma function 

(3.70) T{z;p,q)= [[ ^_ , 

defined for \p\, jqj < 1. This function can be traced back to E.W. Barnes in 1904 [18], 
but was given prominence through the recent work of S.N.M. Ruijsenaars [131]. It 
permits the extension of the standard gamma recurrence to 

^{qz;p,q) = 0{z;p)T{z;p,q), 

where 9(z;p) = {z; p) ao{p / z; p) is a normalised theta function. Another funda- 
mental property of the elliptic gamma function is the functional equation 

(3.71) r(z;p,g)= ^ 



T{pqlz;p,q) 

which follows immediately from the definition ( 3.70[ ). 



The elliptic analogue of the beta integral ( 1.3 1 was discovered in 2000 by V.P. Spi- 
ridonov [137] 

(3-^2) I^^T%^^^-r^^ n r(M.P,.), 
Jc i^yz^ \P,q) 27riz [p;p)ao{q;q)ao i<;^,'^<6 

where each \tj.\ < 1, C is the positively oriented unit circle, rir=i — PQ sj^id 

r(te±™;p,g) :=r(te™;p,g)r(tz-";p,<z). 

The p ^ limit is the well-known Rahman [99, 123] integral, which itself is an 
extension of the Askey- Wilson integral [14] . For the reduction of this last integral 



to the beta integral (1.3 1 see [63] 



J.F. van Diejen and Spiridonov [36] have given an n-dimensional generalization 



of (3.721 which may be viewed as an elliptic extension of the Selberg integral. This 
integral, the p ^ limit of which was first obtained by Gustafson [68], takes the 
form 

.3 73^ / TT r(tzf^zf ;p,g) " nLir(t.^.^^P,g) dz, dz„ 

2"n! 



{p;p)^{q;q)^ j-J-^\T{t;p,q) 
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where |t|, . . . , jigl < 1 and t 



2n-2 



n 



pq. van Diejen and Spiridonov 



provided a proof of (3.73 1 along the Unes of the Anderson and Gustafson proofs of 



(1.1 1 and (2.561 respectively. This required an elliptic generalization of the Dixon- 



Anderson integral (2.241 which, initially, was proved making an assumption about 
the vanishing of certain elliptic integrals. A complete proof of the elliptic Dixon- 
Anderson integral was subsequently given by Rains [124] and by Spiridonov [139]. 



The reduction of the elliptic Selberg integral (3.73 1 to the ordinary Selberg inte 



gral is rather cumbersome, requiring several limits, variables changes and special- 
izations of parameters [126]. Fairly straightforward, however, is to see that (3.73) 



provides an elliptic extension of Gustafson's BC„ constant term identity (2.561. 
To see this one first needs to eliminate Ir using 

n'^i tr = pq- This gives 



rise to several elliptic gamma functions of the form T{pqA) which, by (3.711, may 
be replaced by l/r(l/yl). After these elementary manipulations the p ^ limit 
can be carried out, using that r(z;0,g) — l/(z;q)ac- Finally taking — and 



interpreting the resulting integral as a constant term identity yields (2.56). 

Analogous to (2.59), the integrand of ( |3.73 1 can be used to define an inner 
product. Rains [124] has specified a family of abelian functions which are biorthog- 
onal with respect to this inner product, extending the Rahman-Spiridonov the- 
ory [123, 138] of such functions to the multivariable setting, as well as generalising 
the Koornwinder polynomials and Okounkov BC„ interpolation polynomials [117] 
to the elliptic level. (These functions were independently introduced by H. Coskun 
and Gustafson in [31] without the use of elliptic Selberg type integrals.) Rains 
also extended the integrand of (3.73) analogous to the 2F1 extension (2.49) of the 



Selberg integral, and obtained transformation formulas for the resulting elliptic hy- 
pergeometric integrals. By considering the reduction of his theory to the Selberg 
level Rains obtained, for example, [125] 

Jo Jo .I 



l<i<j<n 



n 



r(a +{2n-i- j) + X, + fij) 



n 



]ti + (2n - I - J + 1) + A, + fij) j-J-^ 



r((j + i)7)r(i + (j + 1)7) 
r(i + 7) 



This integral, which generalises the (3 — j case of Kadell's integral (2.46) is origi 



nally due to Kadell [81] and (for 7=1) L.K. Hua [73]. Kadell's integral (2.46) also 



has an elliptic analogue, which has the feature that the Dotsenko-Fateev integral 



(2.31) is a special case [125]. 



There are other integrals is the literature referred to as elliptic Selberg integrals, 
although they do not contain the actual Selberg integral as a limiting case. These 
integrals arise as solutions to the Knizhnik-Zamolodchikov-Bernard (KZB) heat 
equation for (2n + l)-dimensional 3(2 modules 



27riK 



du d^u 



+ n{n+\)p'{\T)u. 



Here u — u{X,t), p{X,t) ~ i)' {X,t) /'d{X,T) with differentiation with respect to A, 
and '&{X,t) — 9i{nX,T) is a theta function [160]. 
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To describe the relevant solutions to the KZB equations let <J? be the elliptic 
master function 

n 

a>(ti,...,t„;T) = ni?(t„T)-2" Y[ E{U-t,,T)\ 
where E{t^ r) is the eUiptic analogue of t 

The solutions considered by G. Felder, L. Stevens and Varchenko [52] are the linear 
combinations 

w«,rn(A,r) = J«:,™(A,r) + (-l)"+ijK,„,(-A,r), 

where 



/2 1 i I ^ 
. . . , i„; r)0„,„, (a+-|^, r) J| ax{U,T) dh ■ ■ ■ dt„. 



0<t„<---<ti<l 



Whenever necessary this integral is understood in the sense of analytic continuation 
from the region where the exponents in ^i/** have positive real part [52], |i| — 
ti + ■ ■ ■ + tn, 0K,m{t, t) is a theta function of degree k and characteristic m 



and axit,T) = 9{X - t,T)/{6{X,T)E{t,T)). 

In several instances Felder, Stevens and Varchenko found that the "elliptic Sel- 
berg integrals" ■UK,m(A, r) permit closed form evaluations in terms of theta functions 
and ordinary gamma functions. The simplest case of such an evaluation corresponds 
to 



■U2n+2,n+l(A, t) = (27r) 



i/2 -ni 



i(3r.-l) 
4(n + l) 



X Sn 



2(n + l)' n+l'2(n + l) 



0(A,t)"+i 
l'2fn + l)M-l-V 



where 5„ is the Selberg integral (1.1 1. 

The value distribution of logC(l/2 + it) on the critical line. The final topic 
to be reviewed, following J. P. Keating and N.C. Snaith [89], is a link between the 
Selberg integral in its trigonometric form (1.17), and one of Selberg's theorems 
relating to the Riemann zeta function [135]. The latter gives the value distribution 

of logC(l/2 + it) for large t, asserting that for any rectangle B G C 
(3.74) 



lim — 



logC(l/2 + i<) 
t:T <t<2T, ' J- e B 

i log log T 



To relate (3.74) to (1.17), first note that (jl.lSP can be used to evaluate 

n 

(3.75) (He' 



1 + e' 
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where the average is with respect to the eigenvalue PDF (1.9 1 [15]. Three further 
ingredients are then required: the interpretation of ( |3.75 1 as specifying a distri- 
bution in random matrix theory, an hypothesis linking the value distribution of 
logC(l/2 + it) to the value distribution of logA(z), A(z) := 0"=! ('2^P(i^i) ~ ^) 
being the characteristic polynomial for the random matrices, and the large n form 
of (3.75) deduced from ( 1.18| . Regarding the interpretation, note that 

n n 

(3.76) RelogA(-l) = ^logje'"' +1|, ImlogA(-l) " 



2 ^ 

i=l 



It follows immediately that 



S{s - Relog A(-l))5(i - Imlog A(-l)) 



dsdt 



is equal to (3.75). In other words, the characteristic function for the joint distribu- 
tion of the quantities (3.76) is equal to ( 3.75[ ). 

The hypothesis of Keating and Snaith [89], which extends the Montgomery- 
Odlyzko law linking the Riemann zeros to eigenvalues of large complex Hermitian 
random matrices (see e.g., [88]), asserts that the value distribution of log C( 1/2 -I- it) 
for large t will coincide with the value distribution of log A(z), \z\ = 1, for A{z) the 
characteristic polynomial of matrices from the CUE {(3 = 2 case of (1.9l) for large 
n. Further, the value of n in the CUE is to be related to the value of t in ({1/2 + it) 
by n = log t, which ensures that to leading order the density of eigenvalues and 
zeta function zeros is equal. 

Thus the task at hand is to compute the large-n limit of (3.75) with k i-^ 
k/{^ logn)-*^/^, I H- » Z/(^ log rt)^/^, which for /3 = 2 and with the identifications of the 

previous paragraph corresponds to the scaling of log C(l/2 + it) by (^ log log T) 



1/2 



in ( 3.74[ ). Using (1.18) this limit has been computed in [15] as being equal to 
exp(— (fc^ + P)/2). Hence the Selberg integral evaluation in its form ( |1.18 1 implies 



that the joint distribution of the scaled logarithm of the characteristic polynomial 
is equal to exp(— (s^ -I- t^)/2), giving quantitive agreement between the hypothesis 
of Keating and Snaith and Selberg's theorem (3.74). 



The value distribution of |A(z)| for |2| = 1 is also of relevance to zeta function 
theory [89] . The characteristic function of this quantity does not lead to a tractable 
integral. On the other hand, with p(s) a distribution supported on s > 0, knowledge 
of the Mellin transform (complex moments) 



m{x) 



^^p{s)ds 



as a function in the complex plane gives, via the inverse Mellin transform. 



p{s) = 



1 

2tti 



c+i oo 



'^ni{x) dx. 



With p{s) the distribution of values of |A(— 1)|-^ for the CUE, m(x -I- 1) is equal 
to (3.75) with k = Q, il = 2x and thus from (1.18) is given explicitly in terms of 



gamma functions (for a discussion of computing the corresponding inverse Mellin 
transform, see [130]). It should also be remarked that the value distribution of 
A(±l) for K{z) the characteristic polynomial of a random orthogonal or unitary 
symplectic matrix, chosen with Haar measure, is a special case of ( |1.23[ ), and thus 
similarly is an example of the Selberg integral. Keating and Snaith [90] make use 
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of this fact to provide a quantitative link between the value distribution of families 
of i-functions on the critical line and random matrix theory. 
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